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CLAIMS 

1 . Method of determining the velocity V and 
anellipticity r| parameters for processing seismic traces 
in a common midpoint (CMP) gather including an anelliptic 
NMO correction, comprising: 
5 - a preliminary step to define a plurality of nodes 

(dtn, x 0 ) , the said nodes being indicative of parameters 
dtn and t 0 representing the NMO correction for the 
maximum offset and the zero offset travel time in 
hyperbolic coordinates, the said preliminary step being 
10 followed by 

- for each node (dtn, t 0 ) defined in the preliminary 
step, the following steps: 

- for static NMO correction of traces in the CMP 
gather as a function of the values of the said parameters 

15 dtn and x 0 at the node considered, and 

- for calculating the semblance function 
associated with the said NMO correction for the node 
considered; and 

- for each picked time t 0 , a step including 
20 determination of the maximum semblance node (dtn(t 0 ) , 

T 0 (t 0 ) ) , 

- and a final step to convert the dtn(t 0 ) and T 0 (t 0 ) 
parameters so as to obtain the velocity V(t 0 ) and 
anellepticity r| (t 0 ) laws. 

25 2 * Method according to claim 1, wherein the nodes 

are defined during the preliminary step in an analysis 
volume (dtn, x 0 , t 0 ) determined by minimum and maximum 
values respectively [dtn min , dtn max ] [T 0min , Xomaxl and [t 0m in, 
t 0max ] of the dtn, t 0 and t 0 parameters. 
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3. Method according to claim 2, wherein, during the 
preliminary step, a corridor [dtrw, (t 0 ) , dtn max (t 0 )J, 
[tomin(to), Tomax(to)] for changing dtn and t 0 parameters ia 
delimited inside the analysis volume as a function of 

5 plausible velocity V and anellipticity ti values, the 
nodes (dtn, x 0 ) defined for applying the NMO correction 
being then located along the corridor thus delimited. 

4. Method according to any one of the preceding 
claims, further comprising, for each node (dtn, x 0 ) , a 

10 stacking step of the corrected seismic traces, following 
the semblance function calculation step. 

5. Method according to claim 4, wherein the stacking 
of corrected traces is done using only near offset 
traces . 

15 6 - Method according to any one of claims 4 or 5, 

further comprising for each picked time, and following 
the step for determining the maximum semblance node, a 
step of checking that values dtn and x 0 of the maximum 
semblance node correspond to a stacking extreme value for 

2 0 the same values dtn and t 0 . 

7. Method according to any one of the preceding 
claims, further comprising a step of selecting and 
adjusting the pickings obtained, following the step 
implemented for determining the maximum semblance node 
(dtn(to), x 0 (t 0 )) for each picked time t 0 , before the 
conversion step. 

8. Method according to claim 7, wherein the said 
step of selecting and adjusting the pickings comprises a 
step of only retaining pickings dtn and t 0 for which time 
to the highest semblance pickings is greater than a 
predefined value. 
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9. Method according to claim 8, wherein the said 
step . of selecting and adjusting the pickings also 
comprises a step for adjusting the retained pickings dtn 
and x 0 by parabolic interpolations using values about the 

5 said picked values. 

10. Method according to claim 9, wherein the said 
step of selecting and adjust ingt pickings also comprises 
a step of eliminating retained and adjusted pickings dtn 
and x 0 when it is impossible to calculate the Dix interval 

10 velocities between the picking considered and higher 
semblance pickings. 

11. Method according to any one of claims 1 to 10 , 
wherein the processing applied to seismic traces is an 
NMO correction process implementing a static correction 

15 CORRnmq. 

12. Method according to claim 11, wherein, during 
the preliminary step, the NMO corrections CORRnmo are 
calculated for all nodes (dtn, t 0 ) included in the 
analysis volume and all offsets of processed seismic 

2 0 traces. 

13. Method according to claim 12, wherein the NMO 
correction carried out for each node (dtn, t 0 ) consists of 
applying NMO corrections CORRnmo calculated during the 
preliminary step. 

25 14. Method according to any one of claims 11 to 13, 

wherein for a given (dtn, t 0 ) pair, the static NMO 
correction CORRnmq of a seismic trace with offset x is 
carried out according to the following equation: 

CORR NMO (x)= ~ To + L2 + dtn(dtn +2 2T 0 ) x2 ^ ^.^ ^ 

V x max 

3 0 represents the maximum Offset in the CMP gather. 
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15. Method according to any one of claims 1 to 10, 
wherein the processing applied to seismic traces is a 
PSTM migration using a static NMO correction CORR PSTM . 

16. Method according to claim 15, wherein, during 
the preliminary step, the NMO corrections CORR PSTM are 
calculated for all nodes (dtn and x 0 ) included in the 
analysis volume and all migration offsets inside the 
migration aperture. 

17. Method according to claim 16, wherein the NMO 
correction step carried out for each node (dtn and t 0 ) 
comprises, for each offset class, application of the said 
NMO corrections CORR PSTMl calculated during the 
preliminary step on all midpoints inside the migration 
aperture . 

15 18 • Method according to claim 17, wherein the NMO 

correction step carried out for each node (dtn and x 0 ) 
comprises, for each offset class, the stack of the 
corrected midpoints following application of the said NMO 
corrections C0RR PSTM . 

20 19 • Method according to any one of claims 15 to 18, 

wherein, for a given pair (dtn and x 0 ) , the static NMO 
correction C0RR PST m is carried out according to the 
following equation: 

Corr PSTM (x) = -T 0 + R , dtn(dtn + 2T Q )(x~x + hy + R" dtn(dtn + 2T 0 )(x-x-h) 2 

U ~*L V4 + ^ 

2 5 where: 

- x m represents the coordinates of the midpoints, 

- x - Xm represents the migration- aperture PSTM, 

- h is the half source - receiver offset, 

- xmax is the maximum offset and aperture of the 

3 0 migration. 
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20. Method according to any one of claims 14 or 19, 
wherein, during the final conversion step, the parameters 
dtn (to) and (x0) are converted' to the velocity law V(t0) 
according to the following equation: 



V = 



\dtn{dtn + 2T Q )^- 



21. Method according to one of claims 14 or 19, 
wherein, during the final conversion step, the parameter 
x0 (to) is converted to the anellepticity r|(t0) law 



according to tj « ~ 



10 22. Method according to claims 20 and 21, wherein 

parameter dtn is defined with respect to the velocity V 
and anellepticity r\ according to the following equation: 



dtn = 




(1 + 8tOV 2 



15 23. Method according to claim 21, wherein parameter 

xO is defined according to anellepticity t] according to 
the following equation: ■ 



T 0 = 



1 + 8t| 

24 . Method of characterising a velocity field for 
2 0 processing seismic data using a gather of seismic traces 
at common midpoint, wherein, for each travel time to for 
a zero offset, a set of parameters dtn and xO is defined, 
representing the NMO correction for maximum offset, and 
the zero offset travel time respectively, in hyperbolic 
2 5 coordinates . 



